Screening of a charged particle by multivalent counterions in salty water: Strong 

charge inversion 
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Screening of a macroion such as a charged solid particle, a charged membrane, double helix DNA 
or actin by multivalent counterions is considered. Small colloidal particles, charged micelles, short 
or long polyelectrolytes can play the role of multivalent counterions. Due to strong lateral repulsion 
at the surface of macroion such multivalent counterions form a strongly correlated liquid, with the 
short range order resembling that of a Wigner crystal. These correlations create additional binding 
of multivalent counterions to the macroion surface with binding energy larger than fesT. As a result 
even for a moderate concentration of multivalent counterions in the solution, their total charge at 
the surface of macroion exceeds the bare macroion charge in absolute value. Therefore, the net 
charge of the macroion inverts its sign. In the presence of a high concentration of monovalent salt 
the absolute value of inverted charge can be larger than the bare one. This strong inversion of charge 
can be observed by electrophoresis or by direct counting of multivalent counterions. 

PACS numbers: 87.14.Gg, 87.16.Dg, 87.15.Tt 



I. INTRODUCTION 



Charge inversion is a phenomenon in which a charged 
particle (a macroion) strongly binds so many counteri- 
ons in a water solution that its net charge changes sign. 
As shown below the binding energy of a counterion with 
large charge Z is larger than ksT, so that this net charge 
is easily observable; for instance, it is the net charge that 
determines linear transport properties, such as particle 
drift in a weak field electrophoresis. Charge inversion 
is possible for a variety of macroions, ranging from the 
charged surface of mica or other solids to charged lipid 
membranes, DNA or actin. Multivalent metallic ions, 
small colloidal particles, charged micelles, short or long 
polyelectrolytes can play the role of multivalent counteri- 
ons. Recently, charge inversion has attracted significant 
attentionEJij. 

Charge inversion is of special interest for the delivery 
of genes to the living cell for the purpose of the gene 
therapy. The problem is that both bare DNA and a cell 
surface are negatively charged and repel each other, so 
that DNA does not approach the cell surface. The goal 
is to screen DNA in such a way that the resulting com- 
plex is positively . Multivalent counterions can be used 
for this purpose. The charge inversion depends on the 
surface charge density, so the cell surface charge can still 
be negative when DNA charge is inverted. 

Charge inversion can be also thought of as an over- 
screening. Indeed, the simplest screening atmosphere, 
familiar from linear Debye-Hiickel theory, compensates 
at any finite distance only a part of the macroion charge. 
It can be proven that this property holds also in non- 
linear Poisson-Boltzmann (PB) theory. The statement 
that the net charge preserves sign of the bare charge 
agrees with the common sense. One can think that this 
statement is even raoce universal than results of PB equa- 
tion. It was shownErH, however, that this presumption of 



common sense fails for screening by Z-valent counterions 
(Z-ions) with large Z, such as charged colloidal parti- 
cles, micelles or rigid polyelectrolytes, because there are 
strong repulsive correlations between them when they are 
bound to the surface of a macroion. As a result, Z-ions 
form strongly correlated liquid with properties resem- 
bling a Wigner crystal (WC) at the macroion surface. 
The negative chemical potential of this liquid leads to an 
additional " correlation " attraction of Z-ions to the sur- 
face. This effect is beyond the mean field PB theory, and 
charge inversion is its most spectacular manifestation. 

Let us demonstrate fundamental role of lateral corre- 
lations between Z-ions for a simple model. Imagine a 
hard-core sphere with radius b and with negative charge 
—Q screened by two spherical positive Z-ions with radius 
a. One can see that if Coulomb repulsion between Z-ions 
is much larger than UbT they are situated on opposite 
sides of the negative sphere (Fig. la). 
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FIG. 1. a) A toy model of charge inversion, b) PB approx- 
imation does not lead to charge inversion. 

If Q > Ze/2, each Z-ion is bound because the en- 
ergy required to remove it to infinity QZe/(a + b) — 
Z 2 e 2 /2(a+ b) is positive. Thus, the charge of the whole 
complex Q* = —Q + 2Ze can be positive. For example, 
Q* = 3Ze/2 = 3Q at Q = Ze/2. This example demon- 
strates the possibility of an almost 300% charge inversion. 
It is obviously a result of the correlation between Z-ions 
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which avoid each other and reside on opposite sides of 
the negative charge. On the other hand, the description 
of screening of the central sphere in the PB approxima- 
tion smears the positive charge, as shown on Fig. lb and 
does not lead to the charge inversion. Indeed, in this 
case charge accumulates in spherically symmetric screen- 
ing atmosphere only until the point of neutrality at which 
electric field reverses its sign and attraction is replaced 
by repulsion. 

Weak charge inversion can be also obtained as a triv- 
ial result of Z-ions discreteness without correlations. In- 
deed, discrete Z-ions can over-screen by a fraction of the 
"charge quantum" Ze. For example, if central charge 
— Q = —Ze/2 binds one Z-ion, the net charge of the 
complex is Q* — Ze/2. This charge is, however, three 
times smaller than the charge 3Ze/2 which we obtained 
above for screening of the same charge —Ze/2 by two cor- 
related Z-ions, so that for the same Q and Z correlations 
lead to stronger charge inversion. 

Difference between charge inversion, obtained with and 
without correlations becomes dramatic for a large sphere 
with a macroscopic charge Q S> Ze. In this case, dis- 
creteness by itself can lead to inverted charge limited by 
Ze. On the other hand, it-was predictedu and confirmed 
by numerical simulations^!! that due to correlation be- 
tween Z-ions which leads to their WC-like short range 
order on the surface of the sphere, the net inverted charge 
can reach 



Q* = OM^/QZe', (1) 

i. e. can be much larger than the charge quantum Ze. 
This charge is still smaller than Q because of limitations 
imposed by the very large charging energy of the macro- 
scopic net charge. 

In this paper, we consider systems in which inverted 
charge can be even larger than what Eq. (Q) predicts. 
Specifically, we consider the problem of screening by Z- 
ions in the presence of monovalent salt, such as NaCl, 
in solution. This is a more practical situation than the 
salt-free one considered in Ref. ||,|[ Monovalent salt 
screens long range Coulomb interactions stronger than 
short range lateral correlations between adsorbed Z-ions. 
Therefore, screening diminishes the charging energy of 
the macroion much stronger than the correlation energy 
of Z-ions. As a results, the inverted charge Q* becomes 
larger than that predicted by Eq. (Q) and scales lin- 
early with Q. The amount of charge inversion at strong 
screening is limited only by the fact that the binding en- 
ergy of Z-ions becomes eventually lower than ksT, in 
which case it is no longer meaningful to speak about 
binding or adsorption. Nevertheless, remaining within 
the strong binding regime, we demonstrate on many ex- 
amples throughout this work, that the inverted charge, in 
terms of its absolute value, can be larger than the original 
bare charge, sometimes even by a factor up to 3. We call 
this phenomenon strong or giant charge inversion and its 
prediction and theory are the main results of our paper 



(A brief preliminary version of this paper is given in Ref. 

Since, in the presence of a sufficient concentration of 
salt, the macroion is screened at the distance smaller 
than its size, the macroion can be thought of as an over- 
screened surface, with inverted charge Q* proportional 
to the surface area. In this sense, overall shape of the 
macroion and its surface is irrelevant, at least to a first 
approximation. Therefore, we consider screening of a 
planar macroion surface with a negative surface charge 
density —a by finite concentration, N, of positive Z-ions, 
and concentration ZN of neutralizing monovalent coions, 
and a large concentration Ni of a monovalent salt. Corre- 
spondingly, we assume that all interactions are screened 
with Debye-Huckel screening length r s — (8itIbNi) ' , 
where Ib — e 2 / '{DksT) is the Bjerrum length, e is the 
charge of a proton, D ~ 80 is the dielectric constant of 
water. At small enough r s , the method of a new bound- 
ary condition for the PB equation suggested in Ref. ||,|| 
becomes less convenient and in this paper we develop 
more universal and direct theoretical approach to charge 
inversion problem. 

Our goal is to calculate the two-dimensional concentra- 
tion n of Z-ions at the plane as a function of r s and N. 
In other words, we want to find the net charge density of 
the plane 

a* = -cr + Zen. (2) 

In particular, we are interested in the maximal value of 
the "inversion ratio" , a* /a, which can be reached at large 
enough N. The subtle physical meaning of a* should be 
clearly explained. Indeed, the entire system, macroion 
plus overcharging Z-ions, is, of course, neutralized by the 
monovalent ions. One can ask then, what is the meaning 
of charge inversion? In other words, what is the justifica- 
tion of definition of Eq. (|^) which disregards monovalent 
ions? 

To answer we note that under realistic conditions, ev- 
ery Z-ion, when on the macroion surface, is attached 
to the macroion with energy well in excess of fceT. At 
the same time, monovalent ions, maintaining electroneu- 
trality over the distances of order r s , interact with the 
macroion with energies less than fc^T each. It is this 
very distinction that led us to define the net charge of the 
macroion including adsorbed Z-ions and excluding mono- 
valent ions. Our definition is physically justified, it has 
direct experimental relevance. Indeed, it is conceivable 
that the strongly adsorbed Z-ions can withstand pertur- 
bation caused by the atomic force microscopy (AFM) ex- 
periment, while the neutralizing atmosphere of monova- 
lent ions cannot. Therefore, one can, at least in princi- 
ple, count the adsorbed Z-ions, thus directly measuring 
a*. To give a practical example, when Z-ions are the 
DNA chains, one can realistically measure the distance 
between neighboring DNAs adsorbed on the surface. In 
most cases, similar logic applies to an electrophoresis 
experiment in a weak external electric field such that 
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the current is linear in applied field. Sufficiently weak 
field does not affect the strong (above &_bT) attachment 
of Z-ions to the macroion. In other words, macroion 
coated with bound Z-ions drifts in the field as a single 
body. On the other hand, the surrounding atmosphere 
of monovalent ions, smeared over the distances about r s , 
drifts with respect to the macroion. Presenting linear 
electrophoretic mobility of a macroion as a ratio of ef- 
fective charge to effective friction, we conclude that only 
Z-ions contribute to the former, while monovalent ions 
contribute only to the latter. In particular, and most im- 
portantly, the sign of the effect - in which direction the 
macroion moves, along the field or against the field - is de- 
termined by the net charge a* which, once again, includes 
Z-ions and does not include monovalent ones. Further- 
more, for a macroion with simple (e.g., spherical) shape, 
the absolute value of the net macroion charge can be also 
found using the mobility measurements and the standard 
theory of friction in electrolytesEj. This logic fails only 
for the regime which we call strongly non-linear. In this 
regime, majority of monovalent ions form a bound Gouy- 
Chapman atmosphere of the inverted charge, and, while 
surface charge as counted by AFM remains equal a*, 
the electrophoretic measurement yields universal value 
e/2irlBr s , which is inverted but is smaller than a* . For a 
macroion of the size smaller than r s , its size determines 
the maximum inverted charge. 

Now, as we have formulated major goal of the paper, 
let us describe briefly its structure and main results. In 
Sec. II - IV we consider screening of a charged surface by 
compact Z-ions such as charged colloidal particles, mi- 
celles or short polyelectrolytes, which can be modeled as 
a sphere with radius a. We call such Z-ions "spherical". 
Spherical ions form correlated liquid with properties sim- 
ilar to two-dimensional WC (Fig. 2). 




FIG. 2. Wigner crystal of Z-ions on the background of sur- 
face charge. A hexagonal Wigner-Seitz cell and its simplified 
version as a disk with radius R are shown. 

In Sec. II we begin with screening of the simplest 
macroion which is a thin charged sheet immersed in water 
solution (Fig. 3a). This lets us to postpone the complica- 
tion related to image potential which appears for a more 
realistic macroion which is a thick insulator charged at 
the surface (Fig. 3b). We calculate analytically the de- 
pendence of the inversion ratio, <r*/er, on r s in two limit- 
ing cases r s ^> Rq and r s <C Ro, where R = (na/Ze)~ 1 / 2 
is the radius of a Wigner-Seitz cell at the neutral point 



n = c/Ze (we approximate the hexagon by a disk). We 
find that at r s ^>> Rq 

a* /a = 0M(R o /r s ) = 0.83C 1/2 , (C « 1) (3) 
where ( = Ze/irar 2 = (R /r s ) 2 . At r s < R Q 
a* 2tt( 



a V3 In 2 C 



, (C»i)- 



(4) 



Thus a* I a grows with decreasing r s and can become 
larger than 100%. We also present numerical calculation 
of the full dependence of the inversion ratio on £. 



a) b) 

FIG. 3. Two models of a macroion studied in this paper. 
Z-ions are shown by full circles, a) Thin charged plane im- 
mersed in water. The dashed lines show the position of ef- 
fective capacitor plates related to the screening charges, b) 
The surface of a large macroion. Image charges are shown by 
broken circles. 

In Sec. Ill we discuss effects related to finite size of Z- 
ion. It is well knowrO that monovalent ions can condense 
on the surface of a small and strongly charged spherical 
Z-ion. As a result, instead of the bare charge of Z-ions 
in Eqs. (^) and (|J) one should use the net charge of Z- 
ions, which is substantially smaller. Thus, condensation 
puts a limit for the inversion ratio. The net charge grows 
with the radius a of the Z-ion. Therefore, we study in 
this section the case when r s -C a -C Rq and showed that 
the largest inversion ratio for spherical ions can reach a 
few hundred percent. 

Sec. IV is devoted to more realistic macroions which 
have a thick insulating body with dielectric constant 
much smaller than that of water. In this case each Z- 
ion has an image charge of the same sign and magnitude. 
Image charge repels Z-ion and pushes WC away from the 
surface. In this case charge inversion is studied numer- 
ically in all the range of r s or (. The result turns out 
to be remarkably simple: at C < 100, the inversion ratio 
is twice smaller than for the case of the charged sheet 
immersed in water. A simple interpretation of this result 
will be given in Sec. IV. 

In Sec. V and VI we study adsorption of long rod-like 
Z-ions with negative linear charge bare density —r/o on 
a surface with a positive charge density a. (We changed 
the signs of both surface and Z-ion charges to be closer to 
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the practical case when DNA double helices are adsorbed 
on a positive surface.) Due to the strong lateral repul- 
sion, charged rods tend to be parallel to each other and 
have a short range order of an one-dimensional WC (Fig. 
4) . In the Ref. one can find beautiful atomic force mi- 
croscopy pictures of almost perfect one-dimensional WC 
of DNA double helices on a positive membrane. The 
adsorption of another rigid polyelectrolyte, PDDA, was 
studied in Ref. [l6]. Here we concentrate on the case of 
DNA. 




FIG. 4. Rod-like negative Z-ions such as double helix DNA 
are adsorbed on a positive uniformly charged plane. Strong 
Coulomb repulsion of rods leads to one-dimensional crystal- 
lization with lattice constant A. 

It is well known that for DNA, the bare charge den- 
sity, —rjo is four times larger than the critical density 
— ?7ci— f= — DksT/e of the Onsager-Manning condensa- 
tionEJ. According to the solution of nonlinear PB equa- 
tion, most of the bare charge of an isolated DNA is com- 
pensated by positive monovalent ions residing at its sur- 
face so that the net charge of DNA is equal to — rj c . The 
net charge of DNA adsorbed on a charged surface may 
differ from — rj c due to the repulsion of positive monova- 
lent ions condensed on DNA from the charged surface. 
We, however, show that in the case of strong screen- 
ing, r s <C Aq (A = rjc/p), the potential of the surface 
is so weak that the net charge, —77, of each adsorbed 
DNA is still equal to — r\ c . Simultaneously, at r s <C Aq 
the Debye-Hiickel approximation can be used to describe 
screening of the charged surface by monovalent salt. In 
Sec.V, these simplifications are used to study the case of 
strong screening. We show that the competition between 
the attraction of DNA to the surface and the repulsion of 
the neighbouring DNAs results in the negative net sur- 
face charge density — a* and the charge inversion ratio, 
similar to Eq. (|]): 



— = T-[— t r, {Vc<y/r s > 1) 

a h\(r] c ar s ) 



(5) 



Thus the inversion ratio grows with decreasing r s as in the 
spherical Z-ion case. At small enough r s and a, the inver- 
sion ratio can reach 400%. This is larger than for spheri- 
cal ions because in this case, due to the large persistence 
length of DNA, the correlation energy remains large and 
WC-like short range order is preserved at smaller ar s . 
An expression similar Eq. (|^) has been recently derived 
for the case of polyelectrolyte with small absolute value 
of the linear charge density, 770 <C r] c , and strong screen- 



ing (r s <C A) when screening of both the charged surface 
and the polyelectrolyte can be treated in Debye-Hiickel 
approximation!! The result of Ref. ^ can be obtained if 
we replace the net charge r\ c by the bare charge 770 in Eq. 

In Sec. VI we study the adsorption of DNA rods in 
the case of weak screening by monovalent salt, r s Aq. 
In this case, screening of the overcharged plane by mono- 
valent salt becomes strongly nonlinear, with the Gouy- 
Chapman screening length A = e/iTrlsc*) much smaller 
than r s . Simultaneously, the charge of macroion repels 
monovalent coions so that some of them are released from 
DNA. As a result the absolute value of the net linear 
charge density of a rod, 77, is larger than rj c . We derived 
two nonlinear equations for unknown a* and 77. Their 
solution at r s ^> Aq gives: 



a 



iraa 



■ exp 



'In * In A 
a 2i:a 



V = Vc\ 



I \n(r s /a) 
ln(Ao/2vra) 



(6) 



(7) 



At r s ~ Aq we get 77 ~ rj c , A ~ r s and a* /a ~ 1 so that 
Eq. (0) matches the strong screening result of Eq. (|). 
Since 77 can not be smaller than r) c , the fact that 77 ~ r\ c 
already at r s ~ Aq proves that at r s <C Aq, indeed, 77 ~ r/ c 
In Sec. VII we return to spherical Z-ions and derive 
the system of nonlinear equations which is similar to one 
derived in Sec. VI for rod-like ones. This system lets us 
justify the use of Debye-Hiickel approximation for screen- 
ing of overcharged surface ( Sec. II) at r s smaller than r TO , 
where r m = aexp(i?o/l-65a) is an exponentially large 
length. We show that even at r s 3> r rn nonlinear equa- 
tions lead only to a small correction to the power of r s in 
Eq. (|). 

In Sec. I- VII we assume that the surface charges of 
a macroion are frozen and can not move. In Sec. VIII 
we explore the role of the mobility of these charges. Sur- 
face charge can be mobile, for example, on charged liquid 
membrane where hydrophilic heads can move along the 
surface. If a membrane surface has heads with two differ- 
ent charges, for example, and -e, the negative ones can 
replace the neutral ones near the positive Z-ion, thus ac- 
cumulating around it and binding it stronger to the sur- 
face. We show that this effect enhances charge inversion 
substantially. We conclude in Sec. IX. 



II. SCREENING OF CHARGED SHEET BY 
SPHERICAL Z-IONS 

Assume that a plane with the charge density —a is im- 
mersed in water (Fig. 3a) and is covered by Z-ions with 
two-dimensional concentration n. Integrating out all the 
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monovalent ion degrees of freedom, or, equivalently, con- 
sidering all interactions screened at the distance r s , we 
can write down the free energy per unit area in the form 



T = na 2 r s /D - 2nar s Zen/D + F zz + F, 



id • 



(8) 



where the four terms are responsible, respectively, for the 
self interaction of the charged plane, for the interaction 
between Z-ions and the plane, for pair interactions be- 
tween Z-ions and for the entropy of ideal two-dimensional 
gas of Z-ions. Using Eq. (||) one can rewrite Eq. (@) as 



T = Tr(a*) 2 r s /D + F OC p, 



(9) 



where Focp = F c + Fid is the free energy of the same 
system of Z-ions residing on a neutralizing background 
with surface charge density —Zen, which is convention- 
ally referred to as one component plasma (OCP), and 



F c = -ir(Zen) 2 r s /D + F zz 



(10) 



is the correlation part of Focp- The transformation from 
Eq. (||) to Eq. (6) can be simply interpreted as the ad- 
dition of uniform charge densities —a* and a* to the 
plane. The first addition makes a neutral OCP on the 
plane. The second addition creates two planar capacitors 
with negative charges on both sides of the plane which 
screen the inverted charge of the plane at the distance 
r s (Fig. 3a). The first term of Eq. (§) is nothing but 
the energy of these two capacitors. There is no cross 
term corresponding to the interactions between the OCP 
and the capacitors because each planar capacitor creates 
a constant potential, ^>(0) = 2-Ka*r s / D, at the neutral 
OCP. 

Using Eq. ( |Io| ) , the electrochemical potential of Z-ions 
at the plane can be written as u — Zetp(0) + u ic i + /i c , 
where \iu and /i c = dF c /dn are the ideal and the corre- 
lation parts of the chemical potential of OCP. In equilib- 
rium, a is equal to the chemical potential, /i;,, of the ideal 
bulk solution, because in the bulk electrostatic potential 
ijj = 0. Using Eq. M), we have: 



2na*r s Ze/D = -u c + (n b - u id ). 



(11) 



As we show below, in most practical cases the correla- 
tion effect is rather strong, so that u c is negative and 
| /Li c | 3> k B T. Furthermore, strong correlations imply that 
short range order of Z-ions on the surface should be sim- 
ilar to that of triangular Wigner crystal (WC) since it 
delivers the lowest energy to OCP. Thus one can substi- 
tute the chemical potential of Wigner crystal, /iwc: f° r 
a c . One can also write the difference of ideal parts of the 
bulk and the surface chemical potentials of Z-ions as 



ti b -l*id = k B T]n(N s /N), 



(12) 



where N s ~ n/a is the bulk concentration of Z-ions at 
the plane. Then Eq. (|ll|) can be rewritten as 



where N Q = N s exp(—\u. wc \/k B T) is the concentration 
of Z-ions in the solution next to the charged plane, 
which jplays the role of boundary condition for N(x) when 
x — > CBB. It is clear that when N > Nq, the net charge 
density a* is positive, i.e. has the sign opposite to the 
bare charge density — a. The concentration Nq is very 
small because \^wc\/k B T S> 1. Therefore, it is easy to 
achieve charge inversion. According to Eq. (|lj) at large 
enough N one can neglect second term of the right side 
of Eq. (|ll]) . This gives for the maximal inverted charge 
density 



D \jJwc\ 
2irr„ Ze 



(14) 



Eq. ( |l4| ) has a very simple meaning: \\iwc\l Ze is the 
" correlation" voltage which charges two above mentioned 
parallel capacitors with " distance between plates" r s and 
total capacitance per unit area D/(2,Trr s ). 

To calculate the correlation voltage \nwc\/Ze, we 
start from the case of weak screening when r s is larger 
than the average distance between Z-ions. In this case, 
screening does not affect thermodynamic properties of 
WC. The energy per Z-ion e(n) of such Coulomb WC at 
T = can be estimated as the energy of a Wigner-Seitz 
cell, because quadrupole-quadrupole interaction between 
neigbouring neutral Wigner-Seitz cells is very small. This 
gives 

e(n) = -(2 - 8/3Ti)Z 2 e 2 /RD ~ -l.l5Z 2 e 2 /RD, (15) 

where R — (7rn) _1 / 2 is the radius of a Wigner-Seitz cell. 
A more accurate calculatior£3 gives slightly higher en- 
ergy: 

e(n) ~ -l.UZ 2 e 2 /RD = -1.96n 1/2 Z 2 e 2 / D. (16) 

One can discuss the role of a finite temperature on 
WC in terms of the inverse dimensionless temperature 
r = Z 2 e 2 / (RDk B T). We are interested in the case of 
large V. For example, at a typical Zen = a = 1.0 e/nm 2 
and at room-temperature, V — 10 for Z = 4. Wigner 
crystal meltstJ at T = 130, so that for F < 130 we deal 
with a strongly correlated liquid. Numerical calculations, 
however, confirm that at T > 1 thermodynamic prop- 
erties, of strongly correlated liquid are close to that of 
WCcJ. Therefore, for an estimate of \x c we can still write 
F c = ne(n) and use 



n wc 



d [ne(n)} 
dn 



2^2 



-1.65rfc B T = -1.65- 



. Z 2 e 
RD 



(17) 



We see now that (Mwc is negative and |/Uvkc| ^> fegT, 
so that Eq. ( p^ ) is justified. Substituting Eq. JlTl ) into 
Eq. (H), we get a* = 0.83Ze/ \m s R). At r s > R, charge 
density a* -C a, and Zen ~ a, one can replace R by 
Ro = {an/Ze)- 1 / 2 . This gives 



2Trcr*r s Ze/D = k B T\n(N/N ), 



(13) 



a* /a = 0.83C 



1/2 



(C « i), 



(18) 
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where ( = Ze/nar 2 is the dimensionless charge of a Z- 
ion. Thus, at r s >• R or ( <C 1, inverted charge density 
grows with decreasing r s . Extrapolating to r s = 2R 
where screening starts to modify the interaction between 
Z-ions substantially, we obtain a* = OAa. 

Now we switch to the case of strong screening, r s <C R, 
or £ ^> 1. It seems that in this case a* should decrease 
with decreasing r s , because screening reduces the energy 
of WC and leads to its melting. In fact, this is what 
eventually happens. However, there is a range of r s <C R 
where the energy of WC is still large. In this range, as r s 
decreases, the repulsion between Z-ions becomes weaker, 
what in turn makes it easier to pack more of them on the 
plane. Therefore, a* continues to grow with decreasing 
r s . 

Although we can continue to use the capacitor model 
to deal with the problem, this model loses its physical 
transparency when r s <C R, because there is no obvious 
spatial separation between the inverted charge a* and its 
screening atmosphere. Therefore, at r s <C R, we deal di- 
rectly with the original free energy (0). The requirement 
that the chemical potential of Z-ion in the bulk solution 
equals that of Z-ions at the surface now reads 



dF 
dn 



Hid — Hb ! 



where 



F 



2Trar s Zen 
D 



F : 



zz 



(19) 



(20) 



is the interaction part of the total free energy apart 
from the constant self-energy term ira 2 r s /D. According 
to Eq. (|l2]), at large N when 



M6 - Hid = k B T\n(N s /N) < 2irar s Ze/D 



(21) 



we can neglect the difference in the ideal part of the free 
energy of Z-ion at the surface and in the bulk. There- 
fore, the condition of equilibrium (|l^) can be reduced 
to the problem of minimization of the free energy ( |20| ) 
with respect to n. This direct minimization has a very 
simple meaning: new Z-ions are attracted to the surface, 
but n saturates when the increase in the repulsion energy 
between Z-ions compensates this gain. Since this mini- 
mization balances the attraction to the surface with the 



repulsion between Z-ions, the inequality (21) also guar- 
antees that thermal fluctuations of Z-ions around their 
WC positions are small. Therefore, Fzz can be written 
as 



F : 



zz 



E 



(Zef 



(22) 



where the sum is taken over all vectors of WC lattice. At 
r s <C R, one needs to keep only interactions with the 6 
nearest neighbours in Eq. (E2h. This gives 



F = — 



2irar^Zen 



D 



3n^pexp(-A/r s ), (23) 



where A = (2/^/3) 1 ^ 2 n 1 / 2 is the lattice constant of this 
WC. Minimizing this free energy with respect to n one 
gets A ~ r s InC, R ^ (27r/V^) 1/2 r s ln( and 



a 
a 



2n( 



\/3 ln 2 C 



(C»i). 



(24) 



It is clear from Eq. ( p4| ) that at r s <C R, or ( ^> 1 the 
distance R decreases and inverted charge continues to 
grow with decreasing r s . This result could be anticipated 
for the toy model of Fig. la if the Coulomb interaction 
between the spheres is replaced by a strongly screened 
one. Screening obviously affects repulsion between posi- 
tive spheres stronger than their attraction to the negative 
one and, therefore, makes it possible to keep two Z-ions 
even at Q <C Ze. 

Above we studied analytically two extremes, r s 3> R 
and r s <C R. In the case of arbitrary r s we can find a* nu- 
merically. Indeed, minimizing the free energy ( pp| ) with 
the help of Eq. fl22) one gets 



E 



3 + n/r 

8 Tijr s 



,-n/rs 



(25) 



where the sum over all vectors of WC lattice can be eval- 
uated numerically. Using Eq. (Eq) one can find the equi- 
librium concentration n for any given value of £• The 
resulting ratio a* /a is plotted by the solid curve in Fig. 
5. 
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FIG. 5. The ratio a* /a as a function of the dimension- 
less charge £ = Ze/nar^. The solid curve is calculated for a 
charged plane by a numerical solution to Eq. (p^) , the dashed 
curve is the large r 3 limit, Eq. (|l8|). The dotted curve is cal- 
culated for the screening of the surface of the semispace with 
dielectric constant much smaller than 80. In this case image 
charges (Fig. 3b) are taken into account (See Sec. IV). 
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III. CONDENSATION OF MONOVALENT 
COIONS ON Z-ION. ROLE OF FINITE SIZE OF 
Z-ION. 

We are prepared now to address the question of max- 
imal possible charge inversion. How far can a macroion 
be overcharged, and what should one do to achieve that? 
We see below that to answer this questions one should 
take into account the finite size of Z-ions. 

Fig. 5 and Eq. (|24| ) suggest that the ratio a* /a con- 
tinues to grow with growing £. However, the possibilities 
to increase ( are limited along with the assumptions of 
the presented theory. Indeed, there are two ways to in- 
crease ( = Ze/irar 2 , namely to choose a surface with a 
small a or to choose Z-ions with a large Z. The former 
way is restricted because, according to Eq. (|2l|), Z-ion 
remains strongly bound to the charged plane only as long 
as 2irr s crZe/ D 3> k B Ts where 



s = HNJN) 



(26) 



is the entropy loss (in units of k B ) per Z-ion due to its 
adsorption to the surface. This gives for £: 



C « (max = 2Z 2 l B /i 



(27) 



Therefore, the latter way, which is to increase Z, is re- 
ally the most important one. The net charge Z of a Z-ion 
is, however, restricted because at large charge Z of the 
bare counterion monovalent coions of the charged plane 
(which have the sign opposite to Z-ions) condense on the 
Z-ion surfaced. Assuming that Z-ions are spheres of the 
radius a, their net charge, Z, at large Zq can be found 
from the equation 



Ze 2 /aD = k B Tln(Nx >s /Nx) 



(28) 



where Ze 2 /aD is the potential energy of a monovalent 
coion at the external boundary of the condensation at- 
mosphere ("surface") of Z-ion and k B T \n(Ni, s /Ni) is 
the difference between the chemical potentials of mono- 
valent coions in the bulk and at the Z-ion's surface, 
Ni.g ~ Z/a 3 is the concentration of coions at the sur- 
face layer. Eq. (pi) gives 



Using Eq. 



Z =(2a/l B )\n{r s /a). 
and Eq. (p7|), we arrive at 



(r s > a). 



sl B r s 



In- 



(29) 



(30) 



In the theory presented in Sec. II, the radius of Z-ion, a, 
was the smallest length, even smaller than r s . Therefore, 
the largest a we can substitute in Eq. (|3(]) is a = r s . For 
r s = a = 10A and s = 3 we get (max — 4 so that the 
inversion ratio can be as large as 2. 

Since charge inversion grows with increasing a we are 
tempted to explore the case r s < a C flo- To address 



this situation, our theory needs a couple of modifications. 
Specifically, in the first term of Eq. ( |23| ) we must take into 
account the fact that only a part of a Z-ion interacts with 
the surface, namely the segment which is within the dis- 
tance r s from the surface. Therefore, at r s <C a results 
depend on the shape of ions and distribution of charge. 
If the bare charge of Z-ion is uniformly distributed on 
the surface of a spherical ion this adds small factor r s /2a 
to fiwc and the right side of Eq. (p7j) . This gives 



c 



ai.ix — Z 2 l B /sa 



(31) 



One should also take into account that at a > r s 
should be replaced by 

Z = a 2 /r s l B , 



Eq. as. 



(32) 



which follows from the condition that potential at the 
surface of Z-ion Ze 2 /aD — Ze 2 /(a + r s )D is equal to 
k B Thx(N hs /Ni). Substituting Eq. (§|) to Eq. (|l|) we 
find that ( max is larger than that given by Eq. (|3C 
namely 



Cr, 



2a 3 
sl B r s 2 



, (r s < a). 



(33) 



For a = 20A, r s = 10A, I B = 7A and s = 3 we get 
(max — 8 so that the inversion ratio can be as large as 3. 

Let us consider now a special case of the compact Z- 
ion when it is a short rod-like polyelectrolyte of length 
L < R and radius a < r s . Such rods lay at the surface of 
macroion and form strongly correlated liquid reminding 
WC, so that one can still start from Eq. (27). In this 
case, however, Eqs. (^) and (|32|) should be replaced by 
Z ~ Lr/c/e = L/l B . Thus, (max — 2R 2 / sl B r s and can be 
achieved at L ~ R. 

We conclude this section going back to spherical Z-ions 
and relatively weak screening. Until now we used every- 
where the Debye-Hiickel approximation for description 
of screening of surface charge density a* by monovalent 
salt. Now we want to verify its validity. Theory of Sec. 
II requires that the correlation voltage applied to capac- 
itors \nwc\/Ze is smaller than k B T/e. Using Eqs. (|lj) 
and ( fl7j) one can rewrite this condition as Z < R/1.65l B . 
Substituting Z from Eq. (E9h we find that one can use 



linear theory only when r s < r m , where 
r m = a cxp(i?/3.3a). 



(34) 



For a large R/2a, the maximal screening radius of linear 
theory, r m , is exponentially large. Nonlinear theory for 
r s > r m is given in Sec. VII. 



IV. SCREENING OF A THICK INSULATING 
MACROION BY SPHERICAL Z-IONS: ROLE OF 
IMAGES. 

In Sec. II and III we studied a charged plane immersed 
in water so that screening charges are on both sides of the 
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plane (Fig. 3a). In reality charged plane is typically a 
surface of a rather thick membrane whose (organic) ma- 
terial has the dielectric constant D\ much less than that 
of water D\ -CD. It is well known in electrostatics that 
when a charge approaches the interface separating two 
dielectrics, it induces surface charge on interface. The 
potential created by these induced charges can be de- 
scribed as the potential of an image charge sitting on the 
opposite site of the interface (Fig. 3b). At D\ <c D, 
this image charge has the same sign and magnitude as 
the original charge. Due to repulsion from images, Z- 
ions are pushed off the surface to some distance, d. One 
can easily find d in the case of a single Z-ion near the 
charged macroion in the absence of screening (r s = oo). 
The ^-dependent part of the free energy of this system is 



F = AnaZed/D + {Zef/ADd. 



(35) 



Here the first term is the work needed to move Z-ion 
from the surface to the distance d, and the second term 
is the energy of image repulsion. The coefficient 47r (in- 
stead of 2ir) in the first term accounts for the doubling of 
the plane charge due to the image of the plane. The ion 
sits at distance d = do which minimizes the free energy 



of Eq. (35). Solving dF/dd = 0, one gets 




(36) 



In the presence of other counterions on the surface, the 
repulsive force is stronger, therefore one expects that do 
is a little larger than i?o/4. 

To consider the role of all images and finite r s , let us 
start from the free energy per unit area describing the 
system: 



F = — 



AnanZen 



D 

E 



-d/r. 



(Ze) 



n ^ ( Ze ) 2 -n/n 



Ad 2 



(37) 



where, as in Eq. (|22j), the sums are taken over all vec- 
tors of the WC lattice. The four terms in Eq. (^) are 
correspondingly the self energy of the plane, the interac- 
tion between the plane and the Z-ions, the interaction 
between Z-ions (the factor 1/2 accounts for the double 
counting) , and the repulsion between Z-ions and the im- 
age charges (the factor 1/2 accounts for the fact that 
electric field occupies only half of the space). 

At large concentration of Z-ions in the bulk, the dif- 
ference in the ideal parts of the free energy of Z-ion in 
solution and at the surface can be neglected, therefore, 
one can directly minimize the free energy (|37J to find 
the concentration of if-ions, n, at the surface and the 
optimum distance d. The system of equations 



OF 
~dd 



= 



dF 
dn 



= 0, 



(38) 



can be solved numerically and the results are plotted in 
Fig. 5. A remarkable feature of this plot is that, within 
2% accuracy, the ratio a* /a for the image problem is 
equal to a half of the same ratio for the charged plane 
immersed in water (for which there are no images). If 
we try to interpret this result using Eq. (|lj) of the ca- 
pacitor model (Sec. II) we can say that image charges 
do not modify the "correlation" voltage \u W c\/Ze. The 
only substantial difference between two cases is that for 
the thick macroion, instead of charging two capacitors, 
one has to charge only one capacitor (on one side of the 
surface) with capacitance per unit area D/A?:r s 

The fact that image charges do not modify the " corre- 
lation voltage" can be explained quite simply in the case 
of weak screening r s 3> Rq- In this limit, expanding the 
free energy < ffi\ ) to the first order in d/r s , we get 

, N n , 2-Ko 2 d 2n{a*) 2 r s 

F = ne{n) + -Ze4> wc (n, 2d) + — — + ^ 

(39) 

The physical meaning of this equation is quite clear. The 
first two terms are energies of the WC and of its inter- 
action with the image WC (4>wc (n, 2d) is the potential 
of a WC with charge density Zen at the location of an 
image of Z-ion.) The third term is the capacitor energy 
created by the charge of WC and the plane charge. And 
the final term is the usual energy of a capacitor made by 
the WC and the screening atmosphere. 

At a* I a <C 1 minimization of Eq. ( |39| ) with respect 
to d gives the optimum distance do — 0.3ilo, which is a 
little larger than the estimate (|3^). Minimization with 
respect to n gives an equation similar to Eq. (14) 



D \ywc\ 
Anr s Ze 



(40) 



where u\yc differs from the corresponding value in the 
case of immersed plane (Eq. (|l7|)) only by: 



wc 



d_ 

dn 



'-Zeipwc(n,2do] 



(41) 



It is known that ^c(i) decreases exponentially with 
x when x > A/2tt. Since 2do/(A/27r) ~ 1.8, the poten- 
tial ipwc(n, 2d ) cx exp(— 2d 2n/A) and 5u wc /\u wc \ ~ 
(1 - d ^)exp(-2d /(yl/27r)) ~ 0.02. Thus, at r s > 
Rq the chemical potential fj-wc remains practically un- 
changed by image charges. 

In the opposite limit r s <^ Rq one can calculate the ra- 
tio a* /a by direct minimization of the free energy, with- 
out the use of the capacitor model. Keeping only the 
nearest neighbour interactions in Eq. (B7p one finds 



d = r s In 



2ttC 



c 



7< 



\/31n 2 (C 2 /10(d/r s )) 2^3 In 2 £ 



(42) 
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Comparing this result with Eq. ( p4| ) for the case of im- 
mersed plane (no image charges), one gets 



(c* AQ image _ 1 
if* I Ojno image 4 



1 



In 10 



InC 



(43) 



Eq. ( p3| ) shows that in the limit C, — > oo, the ratio 
ct*/(t for the image problem actually approaches 1/4 of 
that for the problem without image. However, due to 
the logarithmic functions, it approaches this limit very 
slowly. Detailed numerical calculations show that even 
at £ = 1000, the ratio pS] ) is still close to 0.5. In prac- 
tice, C can hardly exceed 20, and this ratio is always close 
to 0.5 as Fig. 5 suggested. 

Although at a given £, image charges do not change 
the results qualitatively, they, as we show below, reduce 
the value of Q m ax substantially. As in Sec. Ill, we find 
Cmax from the condition that the bulk electrochemical 
potential of Z-ions can be neglected. When images are 
present, according to Eq. (p7|), one need to replace the 
right hand side of Eq. (pl|) by 27rar s Ze exp(— dp/r s ). Us- 
ing Eq. ((42|), this condition now reads 



A^Z 2 l B /sr s 



(44) 



Using Eq. (^4j) instead of Eq. (|27| and using Eq. fl29|) for 
Z we get Cmax — 5 at r s = a = 10 A and s = 3. Therefore, 
according to the dotted curve of Fig. 5 which was cal- 
culated for the case of image charges, the inversion ratio 
for a thick macroion can be as large as 100%. 



the other hand, the "bare" surface charge of DNA is 
very large, and its corresponding Gouy-Chapman length 
is much smaller than r s . As the result, one needs nonlin- 
ear theory for description of the net charge of DNA. It 
leads to Onsager-Manning conclusion that positive mono- 
valent ions condense on the surface of DNA reducing its 
net charge, — ry, to — rj c = —DksT/e. Far away from 
DNA, the linear theory can be used. When DNA rods 
condense on the plane, we can still use — i] c as the net 
charge density of DNA, because as we will see later, the 
strongly screened potential of plane only weakly affects 
condensation of monovalent ions on DNA. 

Therefore, we can write the free energy per DNA as 



2~K<jr s Lr\ c 
D 



£ 

i — — cx 
i # 



D 



, iA 

K " I T. 



(45) 



where Kq{x) is the modified Besscl function of 0-th or- 
der. The first term of Eq. ( [45| ) describes the interac- 
tion energy of DNA rods with the charged plane, the 
second term describes the interaction between DNA rods 
arranged in one-dimensional WC, the factor 1/2 accounts 
for the double counting of the interactions in the sum. 

Since the function Kq(x) exponentially decays at large 
x, at r s <C A one can keep only the nearest neighbour 
interactions in Eq. (EBl). This gives 



2nar s Lr] c 



-A/r s ) , (46) 



V. LONG CHARGED RODS AS Z-IONS. STRONG 
SCREENING BY MONOVALENT SALT. 

As we mentioned in Introduction the adsorption of long 
rod-like Z-ions such as DNA double helix on an oppo- 
sitely charged surface leads to the strong charge inver- 
sion. In this case, correlations between rods cause paral- 
lel ordering of rods in a strongly correlated nematic liq- 
uid. In other words, in the direction perpendicular to the 
rods we deal with short range order of one-dimensional 
WC (Fig. 4). 

Consider the problem of screening of a positive plane 
with surface charge density, cr, by negative DNA dou- 
ble helices with the net linear charge density — r/ and the 
length L smaller than the DNA persistence length L p so 
that they can be considered straight rods. For simplicity, 
the charged plane is assumed to be thin and immersed in 
water so that we can neglect image charges. Modification 
of the results due to image charges is given later. Here, 
the strong screening case r s < A is considered (A is the 
WC lattice constant). The weak screening case, r s 3> A, 
is the topic of the next section. 

We show below that at r s < A screening radius r s 
is smaller than the Gouy-Chapman length for the bare 
plane. Therefore, one can use Debye-Hiickel formula, 
?/>(0) = 2n(jr s /D, for the potential of the plane. On 



which is similar to Eq. (|23|). To find A, we minimize 
the free energy per unit area, F — nf, with respect to n, 
where n = 1/ LA is the concentration of DNA helices at 
the charged plane. This yields: 



27T0T S 
Vc 



= ^/Ajr s ®ty{-A/r s ). 



(47) 



Calculating the net negative surface charge density, 
— cr* = — rjc/ A + a, wc obtain for the inversion ratio 



ln(?7 c /aT s ) 



{r s < A). 



(48) 



As we see from Eq. ([47]) , the lattice constant A of WC de- 
creases with decreasing r s and charge inversion becomes 
stronger. 

Let us now address the question of the maximal charge 
inversion in the case of screening by DNA. Similar to 
what was done in Sec. Ill, the charge inversion ratio is 
limited by the condition that the electrochemical poten- 
tial of DNA in the bulk solution can be neglected and 
therefore, DNA is strongly bound to the surface. Using 
Eq. ( |46| ) and ([l7]), this condition can be written by an 
equation similar to Eq. ( ^l| ) 

k B Ts <t:2Trar s LT] c /D or rj c /<jr s < 2ttL/sI b , (49) 



9 



where s = ln(N s ^rjNA/Nr>NA) is the entropy loss (in 
units of ks) per DNA due to its adsorption to the sur- 
face. Ns^dna and Noma are correspondingly the three- 
dimensional concentration of DNA at the charged surface 
and in the bulk. Inequality ([l9]) also guarantees that 
WC-like short range order of DNA helices is preserved. 
To show this, let us assume that the left and right nearest 
neighbour rods at the surface are parallel to each other 
and discuss the amplitude of the thermal fluctuations of 
the central DNA along the axis x perpendicular to DNA 
direction (in the limit r s <C A, we need to deal only with 
two nearest neighbours of the central DNA). At x = 0, 
the free energy of the rod is given by Eq. (fl6|). At x ^ 
the free energy of the central DNA is 



/(*) 



1-nar s r\ c L 2Lr\\ 



D 



+ 



D 



■ cosh I — 
2A \ n 



-A/r, 



(50) 



Using Eqs. (|50j) and (|47|), we find the average ampli- 
tude, xq, of the fluctuations of x from the condition 
f(x ) - /(0) ~ k B T. This gives x ^ r s ln(Ae/2nar s 2 L). 
The inequality (p)|) then gives: 



. A 

x < r s In — 

r s 



r s In In — <C A 



r s In . 

crn 



(51) 




FIG. 6. The ratio a* /a as a function of r\a/or a . The solid 
curve is calculated for a charged plane by numerical solution 
to Eq. The dotted curve is calculated for the screening 

of the surface of the semi-space with dielectric constant much 
smaller than 80. In this case image charges are taken into 
account. 

One can numerically minimize the free energy ( flo] ) at 
all r s < A to find a* /a. The result is plotted by the solid 
curve in Fig. 6. 

Let us now move to the more realistic case of a thick 
macroion, so that repulsion from image charges must be 
taken into consideration. As in the spherical Z-ion case, 
image charges push the WC off the surface to some dis- 
tance d. The free energy per DNA rod can be written 
as 



4nar s L V _ d/r3 1 ^ ^M K J^ 
J D 2 ^ D I r s 

i — — OG v 
i ^ 



- y 

2 t— 1 D 

1 — — OQ 



(52) 



Thus, DNA helices preserve WC-like short range order 
when the condition ( ^9| ) is met. 

This condition obviously puts only a weak restriction 
on maximum value of a* /a. At L = L p = 50 nm and 
s = 3, the parameter i] c /ar s can be as large as 75 and, ac- 
cording to Eq. (E^ the ratio a* /a can reach 15. There- 
fore, we can call this phenomenon strong charge inver- 
sion. 

This limit can be easily reached at a very small a. On 
the other hand, if we want to reach it making r s very 
small we have to modify this theory for the case when r s 
is smaller than the radius of DNA. In a way, this is sim- 
ilar to what was done in Sec. 3 for spherical Z-ions. At 
r s <C a one replaces the net charge of DNA, n c by rj c a/r s 
and adds small factor (r s /7r 2 a) 1 / 2 to the first term of Eq. 
([46|). This modification changes only logarithmic term 
of Eq. ( J48| ) and does not change our conclusion about 
strong charge inversion. 



where the three terms on the right hand side are cor- 
respondingly the interaction between the plane and the 
DNA, between the different DNAs and between the 
DNAs and their images. 

The equilibrium distance do and A can be obtained by 
minimizing the free energy per unit area F — nf with 
respect to d and n = 1 /LA: 



dF _ dF 
dd dn 



0, 



(53) 



This system of equations is solved numerically. The re- 
sult for a* /a is plotted by the dotted curve in Fig. 6. It is 
clear that in the case of DNA, at a given value of rj c /ar s , 
image charges play even smaller role than for spherical 
Z-ions. The ratio a* /a in the case of a thick macroion is 
close to 70% of a* /a for the charged plane immersed in 
water, instead of 50% as in Fig. 5 for spherical Z-ions. 

However, like in the case of spherical Z-ions, image 
charges modify the maximal possible value of rj c /<jr s sig- 
nificantly. When images are present, according to Eq. 
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(|52|), one need to replace in Eq. ( f49| ) 2irar s Lr] c / D by 
27rar s Lr] c exp(— do/r s ). Therefore, the condition that the 
bulk ideal chemical potential can be neglected and, there- 
fore, DNA is strongly bound at the surface has the form 



ksTs <C 2irar s Lr] c exp(—do/r s ) 



(54) 



Similarly to what was done above for the problem of 
charged plane immersed in water one can show that Eq. 
( |54| ) guarantees, also, WC-like short range order of DNA 
helices. In the limit r\ c crr s , keeping only the nearest 
neighbour interactions in the free energy (52) and min- 
imizing with respect to d one gets do — r s h\{rj c / Aar s ). 
Substituting do into Eq. (M) we arrive at the final form 
for the condition of Eq. (EM: 



r/ c /crr s < \J&nL/sl B , (L < L p ). 



(55) 



It is clear that the maximal r\ c j ar s and maximal inversion 
ratio grow with L. For L = L p = 50 nm and s = 3, the 
maximal r] c j ar s — 25. Therefore, according to the dotted 
curve in Fig. 6, the inversion ratio for a thick macroion 
a* I a can reach 4. Such inversion can still be considered 
as strong. 

Until now we talked about relatively short DNA, L < 
L p , which can be considered as a rod. For DNA double- 
helices of a larger length (L 3> L p ) the maximum in- 
version ratio saturates at the value obtained above at 
L = L. p . This happens because even a long DNA can not 
be adsorbed at the surface if for L = L p inequality ( |55|) is 
violated. (See the theory of adsorption-desorption phase 
transition, for example, in Ref. pl| ). 

On the other hand, if inequality ( [55| ) holds at L = L pi 
i. e. at r\ c jar s <C ^JiirL p / sIb, the adsorption of a long 
DNA is so strong that DNA lays flat on the charged 
surface. Since repulsion between neighbouring parallel 
DNA is balanced with attraction to the surface, interac- 
tions between parallel DNA helices are so strong that the 
same inequality guarantees WC-like short range order at 
the length scale L p , even though DNA length is much 
larger than L p . One can verify this statement studying 
lateral fluctuations of a DNA segment with length L p 
similarly to the calculation presented above for the prob- 
lem of charged plane immersed in water (See Eqs. ( |50| ) 
and (prl)). Thus, our theory and the plots of Fig. 5 are 
applicable for a long DNA and, therefore, for any flexible 
polyelectrolyte. 

To conclude this section, we would like to mention an- 
other charge inversion problem similar to the problem 
we considered here. Giant charge inversion can be also 
achieved if a single very long DNA double helix screens 
a long and wide positively charged cylinder with radius 
greater or about the double helix DNA persistence length 
(Fig. 7). 




FIG. 7. A long charged worm-like rods spirals around an 
oppositely charged cylinder to screen it. Locally, the picture 
resembles that of an one-dimensional WC 

In this case, an DNA double helix spirals around the 
cylinder. Neighbouring turns repel each other so that 
DNA forms an almost perfect coil which locally resem- 
bles one-dimensional WC. As a result, the cylinder charge 
inverts its sign: density of DNA charge per unit length of 
the cylinder becomes larger than the bare linear charge 
density of the cylinder. At small r s this charge inver- 
sion can be as strong as we discussed above. If cylin- 
der diameter is smaller than DNA persistent length one 
should add elastic energy to the minimization problem. 
This, of course, will make charge inversion weaker than 
for wider cylinders, but still it can be quite large. We 
leave open the possibility to speculate on the relevance 
of these model systems to the fact that DNA overcharges 
a nucleosome by about 15%u. 

A similar problem of wrapping of a weakly charged 
polyelectrolyte around oppositely charged sphere was re- 
cently studied in the Debye-Huckel approximation in 
Ref. |7J. A strong charge inversion was found in this case 
as well. Charge inversion for a charged sphere screened 
by an oppositely charged flexibla.polyelectrolyte was pre- 
viously observed in experimentu and numerical simula- 
tionsa. 



VI. LONG CHARGED RODS AS Z-IONS. WEAK 
SCREENING BY MONOVALENT SALT 



In this section, we consider screening of a positively 
charged plane by DNA rods in the case of weak screening, 
when r s 3> A. We saw in Sec. II that when the screening 
radius is larger than the lattice constant of WC, the ca- 
pacitor model provides a transparent description of the 
charge inversion. Here we adopt this model, too. How- 
ever, we find out that in the case of rods, the inversion 
charge a* is so large that its screening by monovalent salt 
is nonlinear. In other words, at r s 3> A, the capacitors 
described in Sec. II becomes nonlinear. Correspondingly 
in this case one has to use the solution of the nonlinear 
PB equation for the plane potential: 



V>(0) ~ -{2k B T/De) ln(r s /A). 



(56) 
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where A = e/ircr*lB is the Gouy-Chapman length. It is 
shown below that A <C A <C r s so that the use of Eq. 
( |56| ) is justified. 

The weak screening of the plane potential has also an- 
other important consequence. The net charge density of 
DNA, —77, ceases to be equal to to the Onsager-Manning 
critical density —r/ c . The charge of the plane forces DNA 
to release some of monovalent coions condensed on it, so 
that 77 becomes larger than r] c . Thus, in this case, we have 
to deal with a nonlinear problem with two unknowns, 77 
and a*. 

One can find these unknowns from the two following 
physical conditions of equilibrium. The first one requires 
that the chemical potential of positive monovalent ions 
(coions) in the bulk of solution is equal to the chemical 
potential of coions condensed on the surface of DNA rods 
which, in turn, are adsorbed on the plane. The second 
condition requires that the chemical potentials for DNA 
rods in the bulk solution and DNA rods of the surface 
WC are equal. Let us write the first condition as 

k B Tbx = -eV>(0) + ^ \n(A/2va), (57) 
iVi u 

where N\ and N\ tS are the concentrations of monova- 
lent coion in the bulk and at the DNA surface respec- 
tively. The left-hand side of Eq. ( |57| ) is the entropy 
loss and the right-hand side is the potential energy gain 
when monovalent salt condenses on the DNA surface (the 
potential at the surface of DNA is the sum of ip(0), of 
the nonlinear plane capacitor made and the potential of 
the DNA charged cylinder with radius a and the linear 
charge density —77, screened at the distance A/2?r, by 
neigbouring DNA). Far from the charged plane, DNA 
net charge regains its value ~rj c , the condition of equi- 
librium of condensed monovalent coions on isolated DNA 
rod with those in the bulk can be written in a way similar 
to Eq. (|57 



kfiTlii 



2erj c r s 
— — In — 
D a 



(58) 



Excludingln(iVi )S /JVi) from Eqs. @ and |>|) and us- 
ing Eq. (|56|) we can write the first equation for A (which 
represents a*) and 77 as 



7? c ln- 



i] c In -J- + 77 In 



2na 



(59) 



The equality of the chemical potential of DNA in the 
bulk and of DNA condensed on the plane can be written 
in the form similar to Eq. (nj) 



Lijip(0) = \fj, W c\ 



Lr, 



Ni 



Lrj 2 A Lrjl 



D 



In - - 

a 



D 



In* 

a 



(60) 



As in Eq. (JTJ), we see that a "correlation voltage", 
\l^wc\/ Lrj, charges two capacitors consisting of the over- 
charged plane and its screening atmosphere to a finite 



voltage ■0(0) ■ The new second and third terms on the 
right hand side are due to the change in the net charge 
of DNA, when it condenses on the plane. Specifically, 
the second term is the gain in the entropy of monovalent 
salt released and the third term is the loss in the self 
energy of DNA when its net charge changes from — r/ c in 
the bulk solution to —77 at the plane surface. Here A is 
the screening length near the plane surface. (This can 
be seen from the fact that the three-dimensional concen- 
tration of monovalent salt at the surface is of the order 
Ni a ~ <J* /2eA and the corresponding screening length 

r s ,'sur f = (^N lt J B )- 1/2 ~ (2Xe/wa*l B ) 1/2 ^ A.) 

A formal derivation of Eq. (J60J) is given in the end of 
this section. 

The free energy per DNA of the one-dimensional WC 
of DNA rods at the surface can be written similarly to 
Eq. (Eq) with the screening length r s replaced by A, 



/=- 



2tt(t7/A)A 
D 



1 _ 277 (%A 



i— — 00, i^O 



Lr] 2 2ttA 
~D U ~A 



(61) 



This result can be interpreted as the interaction of DNA 
with its Wigner-Seitz cell (a stripe with length L, width 
A and charge density rj/A). 

The chemical potential [iwc can be easily calculated: 



dlnf] Lr) 2 , 2ttA 



(62) 



where n = l/LA is the concentration of DNA at the 
charged surface. 

Substituting Eqs. ©, ©, and © into Eq. @, 
we arrive at the second equation for 77 and A 

27777, In ^ = -772 In - - ?7 2 In — + 277,77 In -• (63) 
A a Zira 



a 



Solving Eqs. (]59j) and (|6J) together with A = ri/(a + <7*), 
we get 



V - Vc\ 



1 \n(r s /a) 
ln(A /27ra)' 



In - ~ Jin -In 

a V a Zira 



(64) 



(65) 



where Aq = i] c /a. 

Eq. ( |65| ) shows that the theory is self consistent: when 
r s ^> A , one has r s S> A 3> A . This justifies the use of 
nonlinear potential for the plane. Eq. d64) demonstrates 
that 77 3> ?7 C as we anticipated. Eq. (j64), of course, is 
valid only if 77 < 770 , where 770 is bare linear charge density 
of DNA. 
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The ratio a* /a can now be easily calculated by substi- 
tuting A = e/na*lB into Eq. One arrives at Eq. (||) 
which shows that the ratio a* /a increases as r s decreases, 
but remains smaller than unity. When r s ~ A$ one finds 
from Eqs. ( |64| ) and j6q ) that r\ ~ n e , A ~ r s ~ v4 , and 
cr*/cr ~ 1, what matches the Eq. ( |48| ) obtained for the 
strong screening limit (r s <C A). 

Let us now present a derivation of Eq. (|o|). To cal- 
culate the free energy of the system we use the standard 
charging procedure described, for example, in Ref. p2|and 
used for DNA in Ref. |23|,^4]- First, let us start by calcu- 
lating the electrostatic free energy of a DNA dissolved 
in solution, which can be written as the work needed to 
charge the DNA up to the bare value 770 per unit length 



f = L 



(66) 



where <p(r]') is the self consistent surface potential of DNA 
when its charge is 77' per unit length. Following Ref. ^J, 
let us divide this charging process in two steps. First, 
the DNA is charged from up to r\ c . In this step, one 
can use for <p(r]') the linear (Debye-Huckel) potential 



2rj K (a/r s ) 
D Kx(a/r s )a/r s 



2r}' 

— ln(r s /a), (r s » a). 



(67) 



In the next step, DNA is charged from r\ c to 77. In this 
step, one has to use nonlinear potential for <fi{r)'). It can 
be written as a sum 



<Krf) 



in- 



De A{rf) 



2r? c . r s 

In — 

D a 



(68) 



where the first term is the contribution of the interval 
2a > r > a of the distances r from the DNA axis. In this 
interval potential can be approximated by that of the 
charged plane with charge density 77' /2ita. It has Gouy- 
Chapman form with the corresponding Gouy-Chapman 
length A(t}') — ar) c /r]' < a. The second term in Eq. ( |6S| ) 
is the contribution of interval 00 > r > 2a, where we deal 
with a cylinder of radius a and linear net charge density 
—77 c . Now, we can calculate the free energy of a DNA 
rod (which is also the chemical potential of DNA in the 
bulk solution, apart from an ideal part): 

f = L <j>tf)dr{ + L / <j>{rf)drf 

JO Jt) c 

Lrfc r s L2rj c rj 



In- 

D a 



D A(?7o) D a 



(69) 



In the Onsager-Manning condensation theory, one can 
think of the last two terms in the above expression as 
the free energy of the condensation layer. 

When DNA rods are adsorbed on the surface of the 
macroion, the role of ?7 C as a border between the linear 



and nonlinear charging regime is played by the net charge 
77. We calculate the total free energy of the system by 
first charging the plane surface to a and DNA to r\ re- 
spectively, then continue charging the DNA from 77 to 
the final value 770 . The first charging process leads to the 
standard contribution 

Lr]ip(0) + fi wc + ^77- In - , 
D a 

to the chemical potential of DNA, where the three terms 
result from, correspondingly, the capacitor energy of the 
screening atmosphere, the correlation energy of DNA and 
the self energy of DNA. The second charging process 
builds up the condensation layer around each DNA and 
gives a contribution 



4>{r))dr] 



277 c 770 
D 



In 



A(7?o) 



2r? c , 



77) In — 

a 



where the nonlinear potential of Eq. (|6^) was used. 

The chemical potential of DNA on the charged surface 
is the sum of the two above contributions: 



L#(0) + fj, wc 



Lrj 2 A L2i] c T] Q a 

In In 

Da D 



AM 



(570-17) In-. 
D a 



(70) 



Equating this expression to the chemical potential of 
DNA in the bulk (Eq. (|69|)) one gets the desired Eq. 
©■ 

So far, we have dealt only with the screening of charged 
surface by DNA double helices which are highly charged 
polyelectrolytes. The situation is simpler if one deals 
with weakly charged polyelectrolytes whose bare charge 
density 770 is much smaller than 77c. In this case, there 
is no condensation on the polyelectrolyte. Therefore 770 
plays the role of the net charge 77c. In the weak screening 
case, r s 3> 770/cr, this brings about small changes in Eq. 
(|60|), which now reads: 

L V0 m = \ f mc\ f^ln- - |ln^) . (71) 
\ D a D a I 



Substituting Eq. (|5f 
for A, we get 



and (|62| ) into Eq. ( |7l| ) , and solving 



A ~ r s exp ( — In 



(r s > 770/cr). 



(72) 



Nonlinear effects are important when A < r s , or when 
the exponent in the above expression becomes less than 
— 1. This gives the minimal r s at which nonlinear effects 
are still important. 



r m = (570/0-) exp(?7c/77o) 



(73) 



As we see, r m is exponentially large at 7^/770 3> 1. 
This makes this weak screening case practically unim- 
portant. At smaller, more realistic value of r s , one can 
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use Debye-Hiickel linear theory to describe the potential 
of the plane. For r s < rfo/cr, this has been done in Ref. ||. 
The result is an expression similar to Eq. (^) with the 
net charge rj c replaced by the bare charge 770 . 



VII. NONLINEAR SCREENING OF A CHARGED 
SURFACE BY SPHERICAL Z-IONS. 

Let us now return to the screening of the charged plane 
by spherical Z-ions in the case when screening by mono- 
valent salt is very weak. Our goal is to understand what 
happens when screening radius is larger than r m (see Eq. 
(p4|)), so that Debye-Hiickel approximation of Sec. II for 
the description of screening of surface charge density a* 
by monovalent salt fails and a nonlinear description is 
necessary. 

The nonlinearity of screening leads to two important 
changes in the theory in Sec. II. First, the monovalent 
coions condense on the surface of the Z-ion and reduce its 
apparent charge. We discussed this condensation in Sec. 
Ill, but used for the net charge of Z-ion the value ob- 
tained for isolated Z-ion in the bulk solution (Eq. (p8|)). 
In this section we call this charge Z c (this quantity plays 
a similar role as rj c in previous section) and save notation 
Z for the net charge of Z-ion absorbed at the charged sur- 
face as a part of the WC. When positive Z-ions condense 
on the negative surface, a fraction of monovalent nega- 
tive ions, condensed on the Z-ions is released. Therefore, 
strictly speaking, Z > Z c . The charge Z c can be found 
from Eq. (E8h, which in the revised notation reads 



e — = k B T]xi — — 
aD N\ 



(74) 



Here, as in Sec. Ill, N l s is the concentration of monova- 
lent negative ions at the external boundary of the conden- 
sation atmosphere of the isolated spherical Z-ion. The 
net charge of a Z-ion in WC, Z, can be found from the 
condition of equilibrium of monovalent negative ions con- 
densed on a Z-ion of the WC and those in the bulk so- 
lution 



Ze " e( 2 -^ + m )=k B Tln N ^ 



aD 



RD 



(75) 



The term in the parentheses is the total potential of the 
plane and other adsorbed Z-ions at the considered Z-ion. 
This potential is the sum of the negative potential of WC 
and the potential due to the positive net charge a* of 
the plane given by E<x@. Excluding k B Tln(N hs /Ni) 



from Eqs. (Q) and (76) we obtain the first equation for 
two unknowns Z and A, which is similar to Eq. (p9): 



Ze 2 2.2 Ze 2 
~~aD RD 



2k B T\n{r s /\) 



aD 



(76) 



To write the second equation for these unknowns we start 
from the condition that the chemical potentials of Z-ion 



at the charged surface and in the bulk solution are equal. 
In the close analogy with Eq. (60) of Sec. VI, we can 
write 



Ze V(0) = \nwc\ + (Z- 



Z c )k B Th\—^- 



(z b 2 -zy 

2aD 



(77) 



The second and third terms on the right-hand side ac- 
count for the fact that monovalent ions are released when 
Z-ions condense on the plane surface (so that their en- 
tropy is gained) and simultaneously the self energy of 
the Z-ion is reduced. Using Eqs. (|l7|), @ and Q) we 
obtain the second equation for Z and A 

1.65(Ze) 2 (Z - Z c )Z 6 e 2 



2k B TZhx 



A 



RD 



aD 



+ 



Solving Eqs. @ and 
Z 
Z 

and 



— ~ 1 + 



2aD 
8|) we get 
0.56 a 



R 



X 



r s exp 



1.65a, N ls 

In 

2R Ni 



(78) 



(79) 



(80) 



Approximating iVi, s as Ni t 
Ni 



Z/a 3 , we get 



Therefore 



A 



In- 



r s exp 



N 



= 21n^ 
a 



1.65a r s 

m — 

R a 



(81) 



Nonlinear effects are important when A < r s , or when 
the exponent in the above expression becomes less than 
— 1. This gives the minimal r s at which nonlinear effects 
are still important 



a exp(i?/1.65a) 



(82) 



which matches the estimate Eq. (Q) obtained from the 
side of the linear regime. 

The ratio a* jo can be easily calculated from Eq. (pQ) 



a n<rl B r s 



n<rl B r s 



exp 



(1) 



1.65a , r s 
In — 

R a 



-1.65a/ R 



cx r s 



-(l + 1.65a/i?) 



(83) 



Once again, this ratio increases as r s decreasest3. Com- 
paring Eq. (83) to Eq. (|l8|), we see that nonlinear effects 
change the exponent in the dependence of a* /a on r s by 
1.65a/i? <§; 1. Taking into account the fact that it is 
important only when r s is greater than an exponentially 
large critical value r m (see Eq. (p2[)), one can conclude 
from this section that, in practical situation, non-linear 
effects in the problem of screening of a charged surface 
by spherical Z-ions are not important. 
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VIII. SCREENING OF A MACROION WITH A 
MOBILE SURFACE CHARGE. 



So far we have assumed that the bare surface charges 
of the macroion are fixed and can not move. For solid 
or glassy surfaces, colloidal particles and even rigid poly- 
electrolytes, such as double helix DNA and actin, this 
approximation seems to work well. On the other hand, 
for charged lipid membranes it can be violated. The 
membrane can have a mixture of neutral and, for ex- 
ample, negatively charged hydrophilic heads. In a liquid 
membrane heads are mobile so that negative ones can 
accumulate near the positive Z-ion and push the neutral 
heads outside (see Fig. 8). Since the background charges 
are now closer to the counterion, one can immediately 
predict that the energy of the WC is lower and charge 
inversion is stronger than that for the case of an uniform 
distribution of negative heads. 




FIG. 8. A Z-ion and its Wigner-Seitz cell with radius R are 
shown. The negative heads are concentrated in the shaded 
area with radius R-. The rest of the Wigner-Seitz cell is 
occupied by the neutral ones. 

To simplify the calculation of the free energy, and gain 
more physical insight in the problem, let us use the same 
transformation as in the beginning of section II, namely 
we simultaneously add uniform planar charge densities 
—a* and a* to the plane. The first addition makes a 
neutral WC on the plane. While the second addition 
creates the two planar capacitors. The free energy can 
be written as the sum of the energy of WC and two ca- 
pacitors, in the same way as Eq. (||). Therefore, a* is 
given by Eq. (|l4| ) . 

We use below the Wigner-Seitz approximation to cal- 
culate \iwc- This approximation gives the energy per ion 
of WC as the energy of one Wigner-Seitz cell and neglects 
the quadrupole-quadrupole interaction between Wigner- 
Seitz cells. It provides 5% accuracy for the energy of the 
standard WC on an uniform immobile background (see 
Eq. (|l6|)). In the case of mobile charges, as one sees from 
Fig. 7, the quadrupole moment of the Wigner-Seitz cell is 
even smaller than that for WC on an uniform background 
with the same average charge density a. Therefore, in the 
case of mobile charge, the accuracy of the Wigner-Seitz 
approximation is even better. 

For simplicity, we assume the Wigner-Seitz cell of a 



counterion is a disk with radius R = (ir/n) 1 / 2 . The neg- 
ative heads concentrate around the counterion and make 
a negative disk with radius i?_ < R and charge density 
— (j— where <r_ = a/nirR 2 . > a. The rest of the cell 
is occupied by neutral heads (Fig. 5). The fraction of 
negative heads f 2 = R 2 _/R? is fixed for each membrane. 
The uniform charge case is recovered when there are no 
neutral heads so that i?_ — R and / = 1. 

Let us consider the weak screening case r s 3> R- Under 
the transformation mentioned above, we add a disk with 
radius R, density —a* to the Wigner-Seitz cell to neu- 
tralize it. Now, the total energy of a Wigner-Seitz cell is 
the sum of the interactions of the Z-ion with two disks of 
radiuses R- and R, the self energy of the two disks and 
the interaction between the disks: 



e(n) = - 
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(84) 



The integrations in the last term are taken over the disks 
with radius i?_ and R respectively. This last term can 
be written as 



(R-)J(R) 



2ira-a*R 3 
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Q(J), (85) 



where Q(f) is a function of / only and can be evalu- 
ated numerically for each value of / (it decreases mono- 
tonically from 8/3 at / = 1 to at / = 0). Using 

3), one gets from Eq. (84): 
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(86) 



The last two terms is the correction to s(n) due to the mo- 
bility of the surface charge. In the uniform limit, / = 1, 
Q(f) = 8/3, these two terms vanish and one gets back the 
usual formula for the energy per ion of WC in Wigner- 
Seitz cell approximation, Eq. jl5|). 

The chemical potential for a counterion in the mobile 
charge case, "wc,m, can be easily calculated as 



n WC,m 



d[ne(n)} (Ze) 2 



dn 



RD 



1 4 2/g(/) \ 

f 37T/ 7T J 

(87) 



Here a is approximated by Zen, because at r s 3> R, the 
ratio a* j a <C 1. 

The ratio between chemical potential [iwc.m for the 
mobile charges and the chemical potential /xv^c for the 
immobile charges has been evaluated numerically as func- 
tion of the fraction f 2 of the negative heads. The result 
is plotted in Fig. 9. 
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FIG. 9. The ratio between the chemical potentials nwc,m 
for the mobile charge case and nwc f° r the immobile charge 
case as a function of the fraction of charged heads f 2 . 

Obviously, as / decreases, ^wc,m grows as expected. 
According to Eq. (|l4|) this means that the inversion ratio 
a* I a grows with decreasing /, too. We do not continue 
the plot in Fig. 9 to very small f 2 because in this case, 
the entropy of negative heads plays important role and 
screening by negative heads of the membtane can be de- 
scribed in Debye-Hiickel approximationcJ. We do not 
consider this regime here. 

Let us now move to the limit of strong screening, 
r s <C R- In this limit, it is more convenient to directly 
minimize the free energy, instead of using the capacitor 
model. Since r s -C R, one needs to keep only the near- 
est neighbour interactions in the free energy. Assuming 
r s <C R- , one can write the free energy per unit area as 

F = -2K*-r.Zen + 3» (ze - -) " e ^tAM , (88) 
V n) A 

where Ze — a/n = cr*/n is the charge of one Wigner- 
Seitz cell. In Eq. (jScf), the first term is the interac- 
tion of Z-ion with the negative background (the disk 
with charge density er_), the second term is the interac- 
tion between neighbouring Wigner-Seitz cells. As usual, 
the quadrupole-quadrupole interaction between Wigner- 
Seitz cells is neglected. 

Minimizing the free energy d88| ) with respect to n, one 
gets A ~ r s ln(/ 2 £) and 

— = _ 2 ![ C , (C>1). (89) 

o V31n 2 (/2 C ) K ' 

where C = Ze/irar 2 . Comparing to Eq. (0), one can see 
that, as in the weak screening case, the inversion ratio 
increases due to the mobility of the surface charge. 

Theory of this section is based on the assumption that 
the charge of Z-ion is so large that it is screened nonlin- 
ear ly by the disk of opposite charge. One can easily gen- 
eralize this calculations to rod-like polyelectrolytes and 
study the role of a similar stripe of positive hydrofilic 
heads attracted by strongly negative DNA. Note that 
the idea of nonlinear concentration of charge in mem- 
branes with two types of heads has been used recently in 
a theory of DNA-cationic lipid complexes^!!. 



IX. CONCLUSION 



We would like to conclude with another general phys- 
ical interpretation of the origin of charge inversion. To 
do so, let us begin with brief discussion of a separate 
physical problem, namely, let us imagine that, instead 
of a macroion, a neutral macroscopic metallic particle is 
suspended in water with Z- and mono-valent ions. In 
this case, each ion creates an image charge of opposite 
sign inside the metal and thus attracts to the metal. 
Obviously, this effect is by a factor Z 2 stronger for Z- 
ions than for monovalent ones. While directly at the 
metal surface, energy of interaction of Z-ion with im- 
age, — (Ze) 2 /4a, is much larger than kgT. Therefore Z- 
ions are strongly bound to the metallic surface, making 
it effectively charged, while monovalent ions are loosely 
correlated with the surface, providing for its screening 
over the distances of the order of r s . We can determine 
the net charge of metallic particle with bound Z-ions us- 
ing the "capacitor model" discussed above. Namely, the 
attraction of the Z-ions to their images plays the role of 
correlation part of the chemical potential /i c and provides 
for the voltage Ze/4a which charges a "capacitor" with 
the width r s between metal surface and the bulk solution. 
This leads to the result that metal surface is charged with 
the net charge density a* = Zen = Ze/(16nar s ). Note 
that metallic particle becomes charged due to interac- 
tions, or correlations, between Z-ions and their images, 
even though the particle itself was neutral in the first 
place. 

Major results of the present paper canJpg now inter- 
preted using a similar language of imagefm. Although 
now we consider a macroion with an insulating body, it 
has some bare charge a on its surface, which leads to 
adsorption of certain amount of Z-ions. The layer of ad- 
sorbed Z-ions plays the role of a metal. Indeed, consider 
bringing a new Z-ion to the macroion surface which has 
already some bound Z-ions. New Z-ion repels nearest 
adsorbed ones, creating a correlation hole for himself. In 
other words, it creates an image with the opposite charge 
behind the surface. Image attracts the Z-ion, thus pro- 
viding for the negative [i c in Eq. ( pT| ) and therefore lead- 
ing to the charge inversion. 

The analogy between the adsorbed layer of Z-ions and 
a metal surface holds only at length scales larger than 
some characteristic length. In WC this latter scale is 
equal to Wigner-Seitz cell radius R. This is why for WC 

[i c {Zef/R (see Eq. (0)). To make |/i c | > k B T, 

small enough radius R is needed. This explains why a 
significant bare charge a is necessary to initiate adsorb- 
tion of Z-ions and to create a metallic layer with images 
which can lead to charge inversion. From formal point of 
view, charge inversion in this case can be characterized 
by the ratio a* /a, as we did throughout the paper, while 
for a neutral metallic particle such ratio is infinite. 
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In this paper, we considered adsorption of rigid Z-ions 
with the shapes of either small spheres or thin rods. The 
concept of effective metallic surface and image based lan- 
guage is perfectly applicable in both cases. It appears 
also applicable to the other problem, not considered in 
this paper, namely, that of adsorption of a flexible poly- 
electrolyte on an oppositely charged dielectric macroion 
surfaceEj. To our mind, this idea was already implicitly 
used in Ref. ||, which assumes that Coulomb self-energy of 
a polyelectrolyte molecule in the adsorbed layer is neg- 
ligible. This means that charge of the polyelectrolyte 
molecule is compensated by the correlation hole, or im- 
age. It is the image charge that attracts a flexible poly- 
electrolyte molecule to the surface. Interestingly, con- 
formations of both the polymer molecule and its image 
change when the molecule approaches the surface. 

A similar role of images and correlations is actually well 
known in the physics of metals. In the Thomas-Fermi ap- 
proximation (which ispsimilar to PB one) the work func- 
tion of a metal is zeroE3 (the work function is an analog 
of fi c ). The finite value of the work function is known to 
result from the exchange and correlation between elec- 
trons. For a leaving electron it can be inteporeted as 
interaction with its image charge in the metalE3. 

We believe that interaction with image or, in other 
words, lateral correlations of Z-ions in the adsorbed layer 
is the only possible reason for a charge inversion exceed- 
ing one Z-ion charge (of course, we mean here purely 
Coulomb systems and do not speak about cases when 
charge inversion is driven by other forces, such as, e.g., 
hydrophobicity) . In the Poisson-Boltzmann approxima- 
tion, when charge is smeared uniformly along the surface, 
no charging of neutral metal or overcharging of charged 
insulating plane is possible. 
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